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Abstract We present simultaneous observations of ^^CO, ^'^CO and C^^O ,7=1—0 emis- 
sion in 11 nearby {cz <1000km s^^) bright infrared galaxies. Both ^^CO and ^'^CO are 
detected in the centers of all galaxies, except for ^"^CO in NGC 3031. We have also de- 
tected C^^O, CS J=2-l, and HCO+ J=l-0 emission in the nuclear regions of M82 and 
M51. These are the first systematical extragalactic detections of ^^CO and its isotopes 
from the PMO 14m telescope. 

We have conducted half-beam spacing mapping of M82 over an area of 4' x 2.5' and major 
axis mapping of NGC 3627, NGC 3628, NGC 4631, and M51. The radial distributions of 
i^CO and ^^CO in NGC 3627, NGC 3628, and M51 can be well fitted by an exponential 
profile. The ^^CO/^'^CO intensity ratio, TZ, decreases monotonically with galactocentric 
radius in all mapped sources. The average TZ in the center and disk of the galaxies are 
9.9±3.0 and 5.6±1.9 respectively, much lower than the peculiar Tl{^2A) found in the 
center of M82. 

The intensity ratios of "CO/C^^O, ^^COI\{CO+ and "CO/CS (either ours or literature 
data) show little variations with galactocentric radius, in sharp contrast with the greatly 
varied TZ. This supports the notion that the observed gradient in TZ could be the results of 
the variations of the physical conditions across the disks. The H2 column density derived 
from C^^O shows that the Galactic standard conversion factor (X-factor) overestimates 
the amount of the molecular gas in M82 by a factor of ~^2.5. These observations suggest 
that the X-factor in active star-forming regions (i.e., nuclear regions) should be lower than 
that in normal star-forming disks, and the gradient in TZ can be used to trace the variations 
of the X -factor 
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1 INTRODUCTION 

It is well known that molecules constitute a significant fraction of the interstellar medium (ISM) and stars 
form from molecular clouds. However, the most abundant molecule, H2, cannot be detected directly in 
typical cold (10-40K) molecular cloud, owing to its lack of a permanent electric dipole moment. The 
next most abundant molecule is ^^CO, which is found to be an e xcellent trace r of H2 due to the low 
excitation temperature (~10K) and critical density (~300cm~'^) (lEvanslll999l) . Generally, the lowest 
transition (J=l— 0) of the rotational lines of ^-^CO and its optically thin isotopic variants (e.g., ^^CO 
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and C^^O) can be used to estimate the column density (A^(H2)) o f molecular clouds in our Galaxy 
dFrerkinget al. Ill982t lYoung & Scovillelll982l IWilson et al. Il2009t) . However, up to now, the column 
density of molecular gas in galaxies is difficult to obtain directly in this way due to the beam dilution 
and the weakness of line emission in CO isotopes. Accurate determination of the H2 column densities 
from CO observations has therefore been a longstanding challenge in external galaxies. 

In our Galaxy, results from independent methods show a tight correlation between the integrated 
^^CO line intensity Ii^qq and N{H2), and the ratios of N(H2) to I\'2r.ci appear to be constant with values 
of (l-5)xlO^°cm-^K -^km~^s across the G alaxy (iBohUn et al. IfigTSTlDickman Ill978l iHunter et al. I 
ll997l:IScoville et al. lll987HYoung & Scoville 11199 ih . This constant value is denoted as the ^^CO-to-Ho 
conversion factor, or the X-factor {X = 7V(H2)//i2co)- It is beyond the scope of this paper to analyze 
in great detail the origin of the empirical ' X-factor' since many studies have shown that the X-factor 
varies with different physical conditions and environm ents, and is influenced by t he CO abundance, 
the gas excitation and the radiative transfer processes (IScoville & Solomonlll974l) . For example, the 
amounts of molecular gas in the metal-poor Magellanic Clouds and M3 1 were cl aimed to be under- 
estimated by up to a factor of ^10 if a sta ndard Galactic X-factor was adopted dMalonev & Black"! 
119881: lAllen&Leau"eux1ll993t Ilsraellll997h . Nevertheless, recent imaging of CO clouds reports simi- 
lar standard X-factor in metal-poor galaxies (iBolatto et al. II2008I) . Conversely, the amounts of molec- 
ular gas in the ultralumino us infrared galaxies (ULIRGs) might be overestimated by factors of ~5 
jPownes & Solomon II 19981 ) using the standard Galactic X -factor 

The IRAS all-sky survey has revealed a large population of infrared bright galaxies with bulk en- 
ergy emission in the fa r-infrared ([Soifer et al. 198 7; Sanders et al. 2003), which is mostly due to dust 
heating from starburst (ISanders & Mirabel 1996 1. Ear ly numerous CO ob servations of infrared galax- 
ies found that these galaxies are rich in molecular gas ( lYoung et al. Ill995h. and there exist a correlation 
between CO and far-infra red luminosity (ISolomon & S age 1988l: lYoung & Scovine 1 fl99ll) . though the 
correlation is non-linear (iGao & Solomon Il2004bl) since Lm/Lco correlates with Lir. Moreover, re- 
cent studies on high-redshift star forrning galaxies further confi rmed the validity of CO-IR luminosity 
correlation (ISolomon & Vanden Bout1l2005UDaddi et al. II2OIOI) . 

Although large quantities of observational studies have aimed at mapping the molecular gas distri- 
bution/kinematics in nearby infrared bri g ht galaxies with single - dish and/or interferome t er telescopes 
dBraine et al. 1 1 19931; lYoung et al. 1 1 19951 ISakamoto et al. Ill999t iNishivama et aTl I2OOII: iHelfer et al. " 



I2OO3I: iLerov et al. I l2009b . these surveys are almost all based on the observations of CO J=l— or 
J=2— 1, and only Umited systema tical studies of CO isotope variants such as ^■^CO, have been pub- 
Hshed so far dPagUone et al. 11200 ll). Owing to the large isotope abundance ratio [i^CO]/["CO]w30-70 
across the Galaxy ( Langer & Penziasll990l) . the opacity of ^'^CO is much lower than that of ^^CO, and 
thus ^'^CO is believed to be mostly optical thin and trace the molecular gas column density adequately 
in most cases. Consequently, the variations in the intensity ratio of ^^CO to ^'^CO as a function of galac- 
tocentric radius, could give a reliable test of whether X-factor varies systematically within galaxies. 
Studies of the Galactic X-factor have revealed large variations in its value towar ds the Galactic nuclear 
region, with an orde r of magnitude increase fro m the center to th e outer disk dSodroskiet al.lll995l; 
iDahmen et al. 1 fT998h . In other galaxies, however. iRickard & Blitz claimed that the value of in- 

tegrated intensity rati o of ^^CO to ^'^CO i n the nucl eus is on average a fac tor of two higher than that in 
the disk. In fact, both lYoung & Sanders 1 dl986l) and fSage & Isbelll dl99ll) suggested that there was no 
clear evidence of a systematic difference in TZ within a galaxy, and they suggested that the variations 
in TZ observed bv lRickard & Blitzl (1 19851) is likely to be caused by pointing errors. Recently, mapping 
surveys of ^^CO and ^'^CO emission towards a large sample of nearby galaxies were carried out by 
iPaghone etand200ll) . finding that the same physical processes (e.g., the molecular gas kinetic temper- 
ature) may both affect TZ and X-factor, moreover, the latter is also expected to decrease from the disks 
to the nuclear regions by factors of 2-5. 

However, any study of the physical properties of molecular gas that involves using line intensity 
ratios might be influenced by measurement errors, owing to the telescopes' different beam sizes, uncer- 
tain beam efficiencies, pointing errors, and calibrations. Recently updated and improved sensitivity and 
system performance with the Purple Mountain Observatory (PMO) 14m telescope at Delingha, China 
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allow us to simultaneously observe extragalactic sources systematically in ^^CO, ^'^CO, and C^^O for 
the first time. Consequently, our simultaneous observations of three CO isotope variants with same tele- 
scope are better suited to obtain the well-calibrated line intensity ratios than observations from different 
telescopes carried out with different tunings at different time. 

In this paper, we present the results of simultaneous observations of ^^CO, ^'^CO, and C^^O along 
the major axes in nearby infrared bright galaxies. The sample selection, observations, and data reduction 
are described in §|2j the observed CO spectra, CO radial distributions, and position-velocity diagrams 
are presented in § [3] These results together with the radial distributions of molecular gas, CO isotopic 
ratio TZ, and possible physical mechanisms that could be responsible to the observed variations in TZ are 
discussed in §21 Finally a summary is presented in §|5] A stability analysis of the PMO 14m telescope 
is presented in Appendix |A] 

2 SAMPLE, OBSERVATIONS AND DATA REDUCTION 

The g alaxies in this study were selected from the Revised Bright Galaxy Sample fRBGS. ISanders et al. I 
|2003[) . We selected 11 galaxies based on the following three criteria: (1) /eo^m > 50 Jy or /loo^m > 
lOOJy. This infrared flux cutoff was chosen to ensure both ^^CO and ^"^CO could be detected with 
reasonable integration time si nce it is well known that the infrared luminosity of galaxies is correlated 
with the CO luminosity (e.g.. lSolomon & Sage L[l988l) . (2) cz < 1000km s^^. This velocity limit was 
chosen due to the limited tuning range of the SIS receiver with the PMO 14m telescope. (3)9h< R.A. 
<15h and Decl.> 0°. For the reason that we can take full advantage of the Galactic dead time to observe 
galaxies in the northern sky. Some general properties of the sample galaxies are summarized in Table [T] 

Our observations were made between 2008 January and 2009 December using the PMO 14m mil- 
limeter telescope at the Delingha. We used the 3mm SIS receiver operated in double-sideband mode, 
which allowed for simultaneous observations of three CO isotope variants, with ^^CO in the upper side- 
band and "CO and C^O in the lower sideband. The Half Power Beam Width (HPBW) is - 60", and 
the main beam efficiency r/mb = 0.67. Typical system temperatures during our runs were about 180-250 
K. The FFT spectrometer was used as back end with a usable bandwidth of 1 GHz and a velocity res- 
olution of 0.16 km s^^ at 115GHz. Observations were done in position switching mode and calibrated 
using the standard chopper wheel method. The absolute pointing uncertainty was estimated to be ~10" 
from continuum observations of planets, and the pointing were checked every two hours by taking spec- 
tra toward the standard source IRCh-10216 throughout our observations. Each galaxy was first observed 
at the center position, and then along its major axis from the center to the outer disk positions sepa- 
rated by half -beam size. Besides the CO observations, we also observed the dense molecular gas tracers 
HCO"*" and CS in the nuclear regions of galaxies. 

The data were reduced using CLASS, which is part of the GILDASQ software package. All scans 
were individually inspected, and those with either strongly distorted baselines or abnormal rms noise 
levels were discarded. Line-free channels which exhibited positive or negative spikes more than 3 a 
above the rms noise were blanked and substituted with values interpolated from the adjacent channels, 
and then a linear baseline was subtracted from the Tine-free' spectrum. After correcting each spectrum 
for main beam efficiency Tymb, the temperature scale of the spectra was converted to the main beam 
temperature Tmb scale from the T^. The spectra were then co-added, weighted by the inverse of their 
rms noises, and the final spectral resolution was smoothed to ~20km s^^. 

Since this is the first time that the PMO 14m millimeter telescope systematically observed galaxies 
other than our own, relatively rather long observing time at each pointing is devoted to accumulate 
relatively adequate integration time. Thus our observations are the first long integrations for the PMO 
14m telescope ever conducted and offer an excellent opportunity to test the stability of the upgraded 
system (see Appendix |A]i. The work present here represents the total of ^ 500 hours observing time in 
CO observations before the telescope system maintenance in the summer of 2009 at the PMO 14m and 
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^ 200 hours after the upgrade, with ^ 40% data discarded owing to the unstable spectrum baseline and 
bad weather conditions. 



3 RESULTS AND ANALYSIS 
3.1 CO Isotopic Spectra 

We detected ^^CO emission from the centers of all 1 1 galaxies observed, of which 10 were also detected 
in i^co simultaneously. Four galaxies (NGC 3627, NGC 3628, NGC 4631 and M51) were observed 
along the major axes, and ^^CO emission was detected at all 42 off-center positions along the major axes, 
while ^'^CO was detected at 27 out of the 42 positions. For M82, the starburst galaxy, ^^CO emission was 
detected at all 47 positions that were mapped in the central 4' x 2.5' region, while ^''CO was tentatively 
detected at 15 positions. C^^O emission was only detected at 13 points close to the nuclear regions in 
M51 and M82. 

Here, we focus on presenting ^^CO and C^^O spectra as compared to ^^CO observed simultaneously 
at the same positions, ^^CO spectra at those ^^CO undetected positions won't be shown her e. Since many 
previo us observations of ^^CO emission in these galaxies are available in literature (e.g.. lYoung et al. I 
119951) . and our CO observations show similar results and comparable spectra. The spectra of both 
^^CO and ^'^CO, as well as Spitzer IRAC 3.6 fim, 8 /im and MIPS 24 /im infrared composite color 
images showing the observed beam positions of the mapped galaxies, are shown in Fig.l. All the three 
CO isotopic spectra at these C^^O detected positions are shown in Fig.|2] 



3.2 Derived Parameters 

The observational results and derived properties are summarized in Table |2] The velocity-integrated 
^^CO (and isotopes) intensity, Ico = / T^hdv, is obtained by integrating Tmh oysr th e veloc- 
ity range with CO line emission feature. Using the standard error formulae in iGao I (119961) (also m 
iMatthews & Gaol 1200 ih . the error in the integrated intensity is 

Al = T„„, At.Fwzi/[/ (1 - AvFwzi/W)]'/^ [K km (1) 

where Tims is the rms noise in the final spectrum in mK, / = Aupvvzi/'^u, where Aupvvzi is the 
Unewidth of the emission feature, dv is the channel bin spacing, and W is the entire velocity coverage of 
the spectrum in units of kilometers per second. For non-detections (only CO isotopes), some spectra are 
further smoothed and found to be marginally detected with signal-to-noise of ^2.5. Otherwise, a 2 cr/ 
upper limits are given based on the estimation by using the expected line width from the detected ^^CO 
hnes at exactly the same position. 

The H2 colum n density and mass surf ace density for galaxies in our sample are derived from the 
empirical relations (INishiyama et al. iboOlh 

iV(H2)[cm-2] = 2 X 1020/i2co[K km s"^], (2) 

and 

E(H2)[M0pc~2] = 3.2/i2co[Kkms~^]cosi, (3) 

where cosi corrects the mass to fa ce-on and a Ga lactic ^■^CO-to-H2 conversion factor 
X=2.0xlO^°cm"2K"ikm~^s is adopted dOame et al. 11200 ih . Obviously, it can be seen from Table |2] 
that the column density in M82 is usually a magnitude higher than that in normal spiral galaxies (the 
rest of the sample). 

Assuming that ^'^CO has the same excitation temperature as ^^CO and the molecular cloud is under 
LTE conditions, then we can calculate the average ^-^CO optical depth from 

rPCO)=ln[l-iii(i!^l- (4) 
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where should be corrected for filling factor, and we can only estimat e an average over all of the 
unresolved clouds in the beam. Using the definition in lWilson et al. I (|2009|) . the H2 column density can 
be derived from the ^"^CO line intensity as 

where the "CO abundance ["CO]/[H2] is 8 x 10^^/60 jFrerking et al. Il982l) and the excitation temper- 
ature, TcK, is taken to be the kinetic temperature of the gas, Tk- Thus, we can estimate Tk by equating 
the column densities derived from both ^^CO and ^"^CO. Both the derived t(^'^CO) and Tk are listed 
in Table [2] Note that LTE assumption is most likely invalid in the central regions of M82. Therefore, 
the extremely low optical depth of ^'^CO (~^0.04) should be treated only as a lower limit, and the result- 
ing kinetic temperature (between 70 and 120 K), which are about 3-4 times higher than that in normal 
galaxies, should be treated as an upper limit. 

Similarly, we estimated the optical depth of C^^O ad opting the same m ethod as equation (4), and 
derived the H2 column density from C^^O intensity using (ISato et al. II 19941) 

where the abundance of [C"0]/[H2] is 1.7x 10"'^ (iFrerkinget al. II 19821) and Tox is adopted from the 
value listed in Table |2] The derived values of t{C^^O) and iV(H2)(Ci®0) for the 13 points of M82 
and M51 detected in C^O are Hsted in Table |3] The average optical depths of C^O in M82 and M51 
are 0.02 and 0.05 respectively, both are about three times lower than that of ^"^CO. Therefore, although 
the optical depth of "CO is moderate (r("CO) --0.3-0.4) in a few galaxies (e.g., NGC 3627 and 
NGC 4631; see Table|2]l, C^^O is always optically thin in all cases here. 

3.3 Distribution of CO Isotopes and Line Ratios 

The distributions of ^^CO and ^'^CO velocity-integrated intensities and their ratio 7?. as a function of 
galacto-central radius are shown in Fig. |3] Note that none of these profiles have been corrected for 
inclination. Both the ^^CO and ^^^CO emission line intensities show an exponential decrease with radius, 
and the ratio TZ decreases from nucleus to the outer disk. 



3.3.1 Radial distributions of ^^CO and ^^CO 

The obvious trends shown in Fig. 3 are that the observed radial distribution of ^■^CO follows that of 
^^CO. The integrated line intensities usually peak in the center of galaxies and fall monotonically to- 
ward larger galactocentric radii. For the barred spiral galaxies NGC 3627 and NGC 4631, however, both 
^■^CO and ^-^CO intensitie s peak at radius of ~ 1 kpc (~^0.5') rather than at the nuclei. The same feature 
has also been presented in lPaglione et al. 1 (l200ll) for NGC 463 1 . Some previous high resolution obser- 
vations of barred galaxies revealed that the molecular gas was concentrated in the c entral regions, with 
secondary peaks at the bar ends due to the bar potential and the viscosity of the gas dRegan et al. II 19991 : 
ISakamoto et 111119991) . 

Similar to the stellar luminosity profiles, the observed ^^CO radial distributions in NGC 3627, 
NGC 3628 and M5 1 could also be well fitted by an exponential fit in R 

liR) = la cxp i-R/Ro), (7) 

where R is the galactocentric distance, Rq is the disk scale length, and Iq is the central integrated 
intensity. The solid curves in Fig. [3] show the least-squares fit to the data excluding the center point 
as the nuclear gas could be a distinct component, yielding Iq = 32.6 K km s^^ and Rq = 3.5 kpc for 
the i^co emission in NGC 3628, Iq = 37.8 K km s^^ and Rq = 2.7 kpc and /q = 3.5 K km s^^ and 
Rq = 3.8 kpc for ^^CO and "CO emission in M51. In NGC 3627, the scale lengths are 2.8 kpc and 
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3.9 kpc for the ^^CO and ^■^CO emission, respectively. However, the power law fit with functional form, 
I = Iq R", is more suitable than exponential fit for the ^^CO distribution in NGC 4631. The fit result 
gives a ~ —1.5 and Iq = 41.1 K km s^^. For the other three exponential fit galaxies, the power law 
could also fit ^^CO distribution almost equally well. Th e exponential ^^CO sca le lengths of 2-4 kpc 
equal to 0.2 R25, which is consistent with the results in lNishivama et alH ( l200lb . The distributions of 
^■^CO and ^•^CO emission along the axis with position angle of 90° in M82 are similar to the other four 
galaxies. The distributions along the major axis are not shown here, since the observations of M82 were 
carried out by mapping an area of 4' x 2.5'. 

3.3.2 The integrated line intensity ratios 

The intensity ratio TZ ranges from 3.3±0.7 to 12.8±4.3 for all 36 points both detected in ^^CO and ^^CO 
emission from normal spiral galaxies, with mean value of 9.9±3.0 and 5.6±1.9 for the central and disk 
regions, respectively. However, the average TZ in M82 is about 2.5 times higher t han that in the nucleus 
of normal spiral galaxies. We use the Equation (2) in I Young & Sanders! (1 1986b to calculate the mean 
value of TZ at each radius. The most prevalent trend is that the TZ drops at larger radius in both the barred 
and non-barred spiral galaxies that we have mapped along the major axis in our sample (Fig. 3). Here 
we note that two points 2.5' away from center) are found to have significantly higher ratios (~^9) in 
M5 1 . These abnormal high values in the disks are tended to be observed once the telescope beams are 
located in the most active star-forming regions along spiral arms. 

The detection of C^^O in M82 and M51 also allow us to estimate the intensity ratio of ^"^CO to 
C^^O, which ranges between 0.9±0.7 and 5.3±2.8 with mean value of 2.9±1.4. The rati os measured 
in the nucleus and disk of M51 are 3.7 and 2.6 respectively, agree well with the results of IVila-Vilaro I 
( l2008i) . which claimed first detect C^^O emission in the center of M51 with ^^CO/C^'^O ratio of 3.6 . 
Also the similar ^' ^ CO/C ^^O ratio values have been found in some starburst galaxies (e.g.. lSage et al. L 
ll99ll:lAaltoe"tanil995h . In addition, our detection of C^^O emission in the off-center regions of M5 1 
represents the first report of such detection for this object by far 

Ongoing observations of dense gas tracers HCO+ and CS, toward the nuclear regions of NGC 4736, 
M82, NGC 3628, and M5 1 by far only yield detection in the latter three galaxies. Here, we also use these 
limited dense gas observations, along with the literature data, to help further analyze the CO isotopic 
results. The intensity ratios of ^"^CO to HCO+ and ^■^CO to CS are found to show little variations 
between starburst and normal galaxies, with average values of 1.3±0.3 and 3.2±1.4, respectively. The 
observed integrated intensity and Une ratios and the Uterature data used are listed in Table |3] 

3.4 Kinematics of CO 

Figure |4] shows the CO position-velocity (P — V) diagrams along the major axes of the galaxies 
NGC 3627, NGC 3628, NGC 4631, and M51, and the P - V diagrams with position angle of 0° 
and 90° in M82 as well. It can be seen in Fig.|4]that P — V diagrams along the major axes tend to show 
a gradual increase of rotation velocity in the inner regions (rigid rotation) and a nearly constant velocity 
in the outer regions (differential rotation). 

For NGC 3627 and M51, the P — V spatial velocity maps of ^'^CO are also shown in Fig. |4] 
Obviously, both ^^CO and ^■^CO share essentially similar kinematics and distribution. At each position 
where observations were made, the distribution of Une intensity and the velocity range for ^^CO and 
^■^CO emission are in good agreement. Accordingly, it could also hint that the line ratios at each position 
derived from our observations are relatively quite reliable. 

Figure |5] shows the CO rotation curve and the variation in line width along the major axis. Using 
the mean velocities, inclination, systemic velocity, and position angle of the major axis that listed in 
Table[Tl and by the assumption that the observed velocity reflect only the systemic motion of the whole 
galaxy and circular rotation, the rotation velocity Vr could be derived via 



Vr = {Vohs - V;ys)/sim cos^ , 



(8) 
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where V^bs is an observed velocity of ^^CO, Vsys is the systemic velocity of the galaxy, i is the incli- 
nation angle, and 6 is the azimuth measured in the disk plane. The peak velocity, Knax, in the rotation 
curve of our sample galaxies , ranges from 120 to 240k m s^^ {Vma.x in spiral galaxies are usually 
between 150 and 300km s' MSparke & Gallagher ll2000l) . 



4 DISCUSSION 

4.1 Radial Distribution of Molecular Gas 



In section 13.31 we show that the surface density of the molecular gas in the galaxies observed in our 
sample can be well fitted both by exponential and power law function. The values indicate that the 
data fitting are at about 85% and 90% confidence level by the exponential and power law function, re- 
spectively. We are therefore unable to distinguish clearly between an exponential and a power law radial 
distribution over the regions observe d here with limited sampl ing a nd low resolution . The s imilar con- 
clusion has been poin ted out in both lYoung & Scovillel (Il982l) and lNishivama et al~l (l200ll) . However, 
IScoville et al. I (1 19831) considered that the exp onential distribution w as more suitable than that of power 
law for M5 1 . The region in M5 1 observed bv lScoviUe etan(ll983b is more extended (~ 8.5') than the 
region 7') observed here, and thus the exponential functional form seems to be better for describing 
the distribution of molecular gas in the whole galaxy. In NGC 3627 and M51, the scale l engths Rq of 



the_2fC0 profiles agree well with the optical K-hmd scale lengths of 3.5 and 3.2 kpc (Reganetal. I 
l200l[) . respectively. The results are in line with the finding first noted bv lYoung & ScovilleT i 1982 ) that 



the large-scale radial distribution of the ^^CO gas follows the light of the stellar disk. In fact, some 
recent high resolution observations also reveal that the ^^CO profiles in a m ajority of galaxies f ollow 
the two-component (bulge and disk) stellar brightness distribution quite well (Regan et al. 200lb. Th e 



CO profiles detected in our observations are in good agreement with that in Paglio ne et al. I OOOlh 
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In NGC 3627 and M51, ^'^CO profiles follow similar exponential distributions as ^^CO. Because of 
insufficient data points with significantly high signal-to-noise, however, the detected ^'^CO data in other 
galaxies are too limited to reliably present more useful information on their distributions. 

Comparing the radial distribution of CO integrated intensities in Fig. [3] with the P — V map of CO 
emission intensities in Fig. HI it is found that the intensities are deficient in the center of NGC 3627, 
NGC 4631, and M51. However, the centrally deficient feature in molecular ring emission is apparent in 
NGC 3627 and NGC 4631, but not in M51. So combining the variations in line width shown in Fig.|5] 
we believe the decrement of CO intensity in the central region of M5 1 is likely as a result of the dilution 
in velocity with much wider line width in the center than that in the outer regions of the central disk. 
On the contrary, the molecular ring emission features with little variations in line width within the bar 
region in NGC 3627 and NGC 4631, are probably either due to orbital disrupti on at the inner Lindblad 
resonance or the central holes with gas exhausted by star formation in the nuclei. [Sakamoto et al. I (1 19991) 
have modeled the P — V diagrams to analyze gas distributions, and found that the central gas hole is 
easier to find from P —V diagrams than from integrated intensity maps due to the velocity information. 

4.2 Variations in the Intensity Ratio of ^^^q to "CO 

Previous studies on the variations in TZ in external galaxies have pointed out that some mergers, 
LIRGs/ULIRGs, and the central regions of circumnuclear starburst galaxies tend to have higher values 
of TZ than that in nor mal spiral galaxies (lAalto et al. Ill99l[ fl995; Ca soli et al. 1ll992LlGreveetdni2009l: 
ISage&Isbeir]ll99lh . of which is similar to that in giant molecular clouds of our Galaxy. Besides the 
enhancement of ^^C caused by nucleosynthesis from type-II SN of massive stars, the deficiency of ^^CO 
caused by isotope-selective photodissociation, and the different distrib utions of ^^CO and ^^CO were 
dispu t ed to be alternative reasons for the high T Z value for a long time jAalto et al. II 19951 ICasoli et al. I 
I1992I [Sage et al. Ill99ll ll^niguchi et al. 19991). However, the single-component model calculation of 



non-LTE CO excitation in iPaglione et al. I (1200 Ih suggested that the variations of kinetic temperature. 
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cloud column and volume density, as well as starburst superwinds might all contribute to the explana- 
tions for the variation in TZ. We here explore the possible causes of the observed variations in TZ. 

4.2.1 Possible causes of variations in TZ 

Our results show that TZ varies not only in the nuclei of various types of galaxies, but also within 
galaxies between nuclear regions and outer disks. The Galactic ^^ C/^'^C abundance ratio w as found to 
range from ~'30 at 5kpc in the Galactic center to ~^70 at 12kpc dLanger & Penzias" I I1990I) . However, 
this isotopic gradient is opposite to the gradient in the molecular abundance ratio TZ (Fig. 3). Therefore, 
the enhancement of ^^C in starburst regions is unlikely to be an appropriate explanation for the measured 
high TZ. 

Some authors argue that the selective dissociation of ^•^CO caused by ultraviolet (UV) field in mas- 
sive star formation regions can stimulate the rat io of TZ, since the rarer isotope is less shielded from 
photodissociation jvan Dishoeck & Black |[l988l) . Consequently, C^^O should be even more dissoci- 
ated by UV photons than ^"^CO due to its lower abundance, and a positive correlation between TZ and 
^■^CO/C^^O intensity ratio would expect to exist if this is available. However, Fig.6b shows a very 
marginal anti-correlation between TZ and ^'^CO/C^^O with a correlation coefficient i?=-0.34. Therefore, 
contrary to the expection, our results reveal a weak anti-correlation, and the deficiency of ^•^CO caused 
by isotope-selective photodissociation could be ruled out. 

In addition to the high TZ measured in M 82, the integrated J=2-l/J=l-0 line ratio was found to 
range between 1.0 and 1.4 dMaoet al.lbOOOl) . revealing the existence of highly excited molecular gas. 



In the PDR model of iMao et al. 1 1 200C ), it was suggested that the bulk of ^^CO emission arises from 
warm diffuse interclumpy medium whereas ^'^CO emission originate in denser cores. In the warm PDRs, 
the optical depth of ^^CO emission from the envelope gas could decrease to a moderate value of t 
result in the corresponding r(^'^CO) << 1. Moreover, the large concentrations of molecular gas in the 
nuclear starburst with high Tk can be excited by strong UV em ission from young massive stars, shock s 
and turbulence caused by supernova remnants, and cosmic rays dPineda et al. l2008tlZhang et al. I2OI0I) . 
The most recent Herschel observations of M82 also suggested that turbulence from stellar winds an d 
supernova may dominate the heating of warm molecular gas in the central region (iPanuzzo et al. II2OI0I) . 
Therefore, it is likely to imply that nonthermal motions produced by the stellar superwinds can broaden 
^^CO lines, thus enhance ^^CO emission as more photons located deeper in the clouds are allowed to 
escape. Furthermore, the significant high value of TZ observed in the spiral arms of M5 1 also demonstrate 
that ^^CO emission can be enhanced in active star-forming regions compared with that in inter-arm 
regions. 

4.2.2 X-factor and dense gas ratio in extreme starburst 

Both theoretical and observational investigat ions have revealed that ^^CO emission is influenced by the 
intrinsic properties of molecular cloud (e.g.. 'Pineda et al. \ 2008 : Goodman et al. L l2009t IShettv et al. L 



I2OIII) . In the magneto-hydrodynamic models of iShettv et al. I (l201 it), ^^CO integrated intensity is found 



to be not an accurate measure of the amount of molecular gas, even for ^'^CO, which may also not be 
an ideal tracer in some regions. In this case, the much lower opacity C^^O can give much more reliable 
constraints on H2 column density than optically thick ^•^CO isotopes. Comparing the H2 column density 
derived from ^^CO with that from C^^O listed in Table |2] and Table [3] we find that the amount of 
molecular gas estimated by standard Galactic X-factor is consistent with that derived from C^*0 in 
M51, whereas is overestimated in M82 by a factor of 2.5, equal to the ratio of TZ between M82 and 
M5 1 . Consequently, our results confirm that the X-factor adopted in starburst active regions should be 
lower than that in normal star-forming regions, and the gradient in TZ can trace the variations in X-factor 
However, surveys of C^^O in a larger sample are required to confirm the relation between the variations 
in TZ and X-factor found in this study. 

The average ratio of ^'^CO to C^^O (~2.9±1.4) derived from our observations in M51 and M82 
indicates that a portion of ^'^CO emission has a moderate optical depth, since the ratio of ^■^CO to 
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C^®0 should be ^1 if both ^^CO and C^^O lines are optically th i n jPenzias et al. 1119811) . This result 



is in line with the two-type cloud model suggested in I Aalto et al. I (Il995 ). in which a large fraction of 



CO emission might originate from denser gas component. S ome previous surveys in dense molecule 
have provided support fo r the presence of such dense gas (e.g. JSage et al. iri990HNguyen et ai7lll992l : 
iGao & Solomon L l2004al). In a ddit ion, our detec t ion of HCO+ and CS in M82 and M5 1 is consistent 



with that in lSage et al. I ( Il990l) and lNayloret al. I ( |2010|). The ratios of "CO to HCO+ a nd "CO to CS 



in (U)LIRGs NGC 3256, NGC 6240 and Arp 220 (ICasoli et al. ll992l:lGreve et al. l2009b are found to be 



lower than those in the nuclear regions of normal spirals and M82 (Fig.|6ll, which probably indicate that 
the dense gas fraction is higher for these (U)LIRGs, since ^"^CO can be considered as a relative reliable 
tracer of total molecular gas due to its low abundance. This result agrees well with the conclusion in 
iGao & Solomon I (l2004bl) that galaxies with high star formation efficiency tend to have higher dense gas 
fraction. Moreover, the ratio TZ measured in NGC 3256, NGC 6240 and Arp 220 are found to be much 
higher than those in M82 (Fig. |6]l, which is likely to imply that the bulk of ^^CO emission arise from 
warm diffuse gas is enhanced in the extreme starburst. 

Summarizing the above analysis, the systematical gradient in TZ can be explained by the variations 
in the physical conditions of molecular gas. The standard Galactic X-factor used in M82 overestimates 
the amount of molecular gas by a factor of 2.5, and the variations in X-factor can be well traced by 
the gradient in TZ. Nevertheless, additional observations of both ^^CO and ^^CO lines at J > 2 and 
more C^*0 lines are required to better constrain the physical conditions of the molecular gas in external 
galaxies. 



5 SUMMARY 

We observed simultaneously the ^^CO, ^'^CO and C^^O emission lines in 1 1 nearby infrared-brightest 
galaxies, of which four (NGC 3627, NGC 3628, NGC 4631 and M51) were mapped with half-beam 
spacing along the major axes and M82 was fully mapped in an area of 4' x 2.5'. These are the first 
systematic extragalactic observations for the PMO 14m telescope and the main results are summarized 
as follows: 

l.We detected the ^^CO emission towards 99 of the positions observed, with the ^■^CO seen towards 
51 of these. C^^O was detected at 13 positions close to the nuclear regions in M51 and M82, among 
which the off-center positions in M51 were the first C^*0 detection reported here. 

2. In the four galaxies with major axes mapping, the ^■^CO line intensity decrease from center to outer 
disk, similar to that of "CO. In NGC 3627, NGC 3628, and M51, the radial distribution of both "CO 
and ^'^CO can be well fitted by an exponential function, whereas the ^^CO distribution in NGC 4631 
is better fitted by power law. The scale length of ^^CO emission is about 0.2 R25 with a mean value 
of ~^3 kpc. Moreover, the ^^CO scale lengths in NGC 3627 and M51 are in good agreement with the 
optical scale lengths. 

3. The peak velocity of ^^CO rotation curves ranges from 120 to 240km s^^, and the line widths of 
^^CO lines tend to drop with radius from center to outer disk in all mapped galaxies. Of all positions 
observed, the distribution of both line intensity and profiles of ^^CO and ^"^CO have good agreement, as 
expected with simultaneously ^^CO and ^-^CO observations. Thus, a reliable line intensity ratio TZ can 
be obtained. 

4. A decreasing tendency of TZ with radius from center to outer disk is found in mapped galaxies. 
TZ varies from 3.3±0.7 to 24.8±2.5 in positions with both ^^CO and ^■^CO detected. The average TZ are 
9.9±3.0 and 5.6±1.9 in the center and disk regions of normal spiral galaxies, respectively. 

5. The high TZ measured in M82 is likely to be caused by enhanced ^^CO emission from deeper 
cloud with broad ^^CO line produced by the stellar winds and supernova. The low values of ^'^CO/C^^O 
(^2.8±1.2) found in M82 support the suggestion that a considerable fraction of ^'^CO emission origi- 
nates in denser gas component. Comparing the ratios of ^•^CO/HCO+ and ^"^CO/CS in normal galaxies 
with those in U/LIRGs, the lower values found in U/LIRGs agree with the notion that the galaxies with 
high SFEs tend to have higher dense gas fraction. 
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6. Comparing with the H2 column density derived from C O, the standard Galactic X-factor is 
found to overestimate the amount of molecular gas in M82 by a factor of ~2.5. This confirms the 
assertion that a lower X-factor should be adopted in active starburst regions than that in normal star- 
forming disks, and moreover, the gradient in TZ can be used reliably to trace the variations of the X- 
factor. 
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Fig. la Spectra of ^^CO {thick lines) and ^'^CO {thin lines) from the central region of 
M82 with ^'^CO tentatively detected. All spectra are on the Tmb scale and binned to a velocity 
resolution of ~20km s^^ (for some weak ^"^CO emission positions, the spectra are further 
smoothed to ~40km s^^ for display). ^^CO spectra are divided by 20 for comparison pur- 
poses. The offset from the center position is indicated in each box. A linear baseline has been 
subtracted using the line-free portions of each spectrum. M82 was mapped in the 4' x 2.5' 
central region. The top panel shows the observed positions(crosses) and ^^CO contours (con- 
tours begin at 10 K km s~^ and increase by 40 K km s~^ each step) overlaid on infrared 
image taken from Spitzer {8.0 /xm {red), 5.8 ^m {green), 3.6 ^m {blue)) 
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Fig. lb Same as figure la, but for NGC 3627. The spectra were measured along the major 
(PA=n6°) axes of the galactic disk. The left panel shows the half-beam spacing observations 
overlaid on Spitzer infrared image (24.0/im (red),S.Qfj,m (green),3.6fj,m (blue)). The size of 
HPBW (--60") is represented by circle. 
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Fig. Ic Same as figure lb, but for NGC 3628. The spectra were measured along the major 
(PA=103°) axes of the galactic disk. 
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Fig. Id Same as figure lb, but for NGC 4631. The spectra were measured along the major 
(PA=88°) axes of the galactic disk. 
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Fig. le Same as figure lb, but for m51. The spectra were measured along the major {PA=Q°) 
axes of the galactic disk. 



^CO, ^^CO and C^*0 observations in nearby galaxies 



17 



-0.05 




-200 200 400 
Velocity (km/s) 




300 400 500 600 
Velocity (km/s) 



Fig. 2 Spectra of ^^CO, ^'^CO and C^^O obtained from the same positions on the galaxy M82 
and M51. ^^CO spectra are divided by 30 and 15 for display in M82 and M51, while ^^CO 
spectra are divided by 3 for display. M5 1 disk spectra represent the average emission over the 
disk region except the center (0,0). 
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Fig. 3 Radial distributions of ^^CO and ^'^CO integrated intensity and their ratios at each 
position along the major axes (for M82,the distributions along the axis with position angle of 
90 degree are shown). Error bars are 1 a statistical uncertainty based on the measured rms 
noise in each spectrum. Upper limits (2(t) are denoted with downward arrow for the non- 
detection ^^CO emission, the corresponding lower limits of line ratio are denoted with 
upward arrow. The solid line in NGC 3627, NGC 3628 and M5 1 represent the exponential fit 
to the radial distribution of mean intensity of ^^CO and ^■^CO emission, whUe in NGC 463 1 



^CO, ^^CO and C^*0 observations in nearby galaxies 19 




200 400 200 400 

Velocity (km/s) Velocity (km/s) 




600 800 1000 600 800 1000 600 800 1000 

Velocity (km/s) Velocity (km/s) Velocity (km/s) 




Fig. 4 Position-velocity diagram of ^^CO and ^^CO (only for NGC 3627 and M5 1) emission 
along the major axes of galaxies. For M82, the P —V diagrams with position angle of 0° and 
90° are shown. All the spectra were smoothed to have a velocity resolution of ~20km s^^. 
The dashed contour on the ^'^CO panel represents the lowest ^^CO contour for comparison. 
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Fig. 5 Rotation velocity derived from CO emission line using eq.(8) and line width mea- 
sured at each position along the major axis. The different symbols represent the same as in 
Fig.3 
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Fig. 6 The relationship for the normal spiral and starburst galaxies between the integrated 
intensity ratio of I^CO/^CO and a) i^cO/C^^O, b) ^^CO/C^^O, c) i3CO/HCO+, and d) 
i^CO/CS. The black symbol represents the emission of ^^CO, ^^CO, and C^^O in M82 and 
M51 (add some new detections in IC 342 and NGC 6949); HCO+ and CS in M82, NGC 3628 
and M51. Lower limits (la) of ratio are denoted with right pointing arrow for some non- 
detection C^^O emission. The triangle symbol represents the data taken from literatures (see 
Table gji, the red triangles: NGC 3256, NGC 6240, Arp 220; the blue triangles: NGC 253, 
NGC 1808, NGC 2146, NGC 4826 and Circinus. Note that all the ratios have been corrected 
for the different beamsizes. 
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Table 1 Source List and Galaxy 
Properties 



Source 


Alias 


R.A. 


Decl. 


V 


i 


P.A. 


Type 


D26 


D 


d 






(J2000.0) 


(J2000.0) 


(km s"^) 


(deg) 


(deg) 




(arcmin) 


(Mpc) 


(kpc) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


NGC 2903 


UGC 5079 


09 32 10.5 


+21 30 05.0 


556 


61 


17 


SAB(rs)bc, HU 


12.6x6.0 


7.96 


1.9 


NGC 3031 


M 81 


09 55 33.6 


+69 03 56.0 


-34 


58 


157 


SA(s)ab, LINER/Syl.8 


26.9x14.1 


1.00 


0.3 


NGC 3034 


M 82 


09 55 49.6 


+69 40 41.0 


203 


80° 


65° 


10, Sbrst/HU 


11.2x4.3 


2.90 


1.3 


NGC 3521 


UGC 6150 


1 1 05 49.2 


-00 02 15.0 


805 


58 


164 


SAB(rs)bc, HU/LINER 


11.0x5.1 


11.47 


2.6 


NGC 3627 


M66 


11 20 15.3 


+ 12 59 32.0 


111 


63 


176 


SAB(s)b, LINER/Sy2 


9.1x4.2 


10.41 


2.6 


NGC 3628 


UGC 6350 


11 20 16.2 


+ 13 35 22.0 


843 


89" 


103° 


SAb pec sp, HU/LINER 


14.8x3.0 


12.07 


3.1 


NGC 4631 


UGC 7865 


12 42 07.1 


+32 32 33.0 


606 


85° 


88° 


SB(s)d 


15.5x2.7 


8.67 


2.5 


NGC 4736 


M94 


12 50 52.9 


+41 07 15.0 


308 


35 


100 


(R)SA(i)ab, Sy2/LINER 


11.2x9.1 


4.40 


1.4 


NGC 5055 


M63 


13 15 49.5 


+42 01 39.0 


504 


55 


105 


SA(rs)bc, HU/LINER 


12.6x7.2 


7.21 


2.3 


NGC 5194 


M51a 


13 29 53.5 


+47 11 42.0 


463 


20° 





SA(s)bc pec, HU/Sy2.5 


11.2x6.9 


6.62 


2.2 


NGC 5457 


M 101 


14 03 09.0 


+54 21 24.0 


241 


27 


40 


SAB(rs)cd 


28.8x26.9 


3.45 


1.4 



From lYoung et al. Hl995l) . 

Notes: Cols.(l) and (2): Galaxy name. Col. (3) and (4): Adopted tracking center. Units of right ascension are hours, minutes, and 
seconds, and units of declination are degrees, arcminutes, and arcseconds. Col. (5): Heliocentric velocity drawn from the litera- 
ture and NED. Cols. (6) and (7): Inclination (i) and position angle (P.A.) from Heifer et al.2003, except where noted. Col. (8): 
Morphological type and nuclear classification from the NED database. Col. (9): Major- and minor- axis diameters from the NED 
database. Col. (10): Luminosity distance, calculated assuming Hq = 70 km s^^MPc^^ . Col.(l 1): Linear scale of 60" at distance 
D, taken from the NED database. 
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Table 2 Observed and Derived Properties of ^^CO and 
"CO 



Il2r 



Il2r 



iV(H2)^ 



(arcmin) 


(aixmin) 


(Kkms~l) 


(kms-1) 




(Kkms-1) 


(kms~^) 


(kms^^) 






(10^" 


(K) 


0.00 


0.00 


24.89it0.92 


548 


155 


1.95 ±0.58 


591 


151 


12.8±4.3 


0.08 


49.8 


55 


0.00 


0.00 


0.66±0.28 


-38 


127 


< 0.68 






>1.0 




1.3 




0.00 


0.00 


224.61 ±1.51 


159 


138 


10.25±0.82 


126 


123 


21.9±1.9 


0.05 


449.2 


98 


0.00 


0.50 


163.0it 1.62 


180 


162 


1 1.80±1.06 


153 


167 


13.8±1.4 


0.08 


328.8 


66 


0.00 


-0.50 


113.4±1.12 


150 


1 15 


5.59±0.90 


137 


150 


20.3±3.5 


0.05 


230.0 


1 10 


0.50 


0.00 


344.77it2.04 


214 


184 


13.91±1.31 


230 


194 


24.8±2.5 


0.04 


689.5 


1 1 1 


0.50 


0.50 


278.07itl.53 


233 


183 


1 1.37±0.97 


238 


180 


24.5±2.2 


0.04 


556.1 


1 10 


0.50 


-0.50 


1 34.00 zt 1.50 


169 


125 


5.61±!.10 


162 


87 


23.9±4.9 


0.04 


269.8 


1 10 


1.00 


0.00 


228.83itl.66 


242 


180 


1 1.29±1.25 


265 


142 


20.3±2.4 


0.05 


457.7 


90 


1.00 


0.50 


218.00itl.69 


268 


162 


6.38±0.83 


313 


110 


34.2±4.7 


0.03 


435.8 


120 


1.50 


0.50 


92.34it2.04 


255 


164 


6.99±1.01 


264 


215 


13.2±2.2 


0.08 


184.7 


57 


1.50 


-0.50 


48.95it 1.74 


173 


94 


3.64±1.36 


158 


82 


13.4±5.5 


0.08 


97.9 


58 


2.50 


0.00 


48.64it2.44 


219 


147 


1.94±0.96 


158 


46 


25.1±13.7 


0.04 


97.3 


1 13 


-0.50 


0.00 


75.91it 1.88 


144 


1 12 


6.78±2.83 


1 14 


100 


1 1.2±5.0 


0.09 


151.8 


48 


-0.50 


-0.50 


45.75it 1.66 


140 


1 10 


5.75±1.33 


226 


250 


7.9±2.1 


0.13 


91.5 


32 


-1. 00 


0.00 


47.85it 1.96 


135 


104 


2.52± 1.49 


141 


44 


18.9± 12.0 


0.05 


95.7 


84 


-1.50 


0.00 


47.40it3.86 


162 


147 


5.81±2.16 


183 


111 


8.2±3.7 


0.13 


94.8 


33 


0.00 


0.00 


22.00 itO. 37 


770 


231 


2.05±0.66 


770 


210 


10.7±3.7 


0.10 


44.0 


45 


0.00 


0.00 


17.24it0.87 


769 


179 


1.66±0.49 


866 


200 


10.4±3.6 


0.10 


34.5 


44 


0.00 


0.50 


25.26zt0.83 


653 


157 


4.06±0.98 


625 


163 


6.2±1.7 


0.18 


50.5 


24 


0.00 


1.00 


13.88itO.38 


633 


177 


2.28±0.51 


597 


138 


6.1±1.5 


0.18 


27.8 


23 


0.00 


1.50 


8.55it0.52 


554 


58 


1.46±0.48 


580 


75 


5.9±2.3 


0.19 


17.1 


22 


0.00 


2.00 


6.90it0.50 


583 


75 


1.79±0.35 


527 


74 


3.9±1.0 


0.30 


13.8 


13 


0.00 


-0.50 


22.25 zt 1.96 


827 


1 15 


2.87±0.69 


829 


145 


7.7±2.5 


0.14 


44.5 


31 


0.00 


-1.00 


8.96zt0.53 


851 


185 


1.78±0.35 


885 


91 


5.0±1.3 


0.22 


17.9 


18 


0.00 


-1.50 


6.82zt0.31 


872 


56 


1.96±0.41 


877 


90 


3.5±0.9 


0.34 


13.6 


1 1 


0.00 


-2.00 


5.90zt0.42 


886 


66 


1.02±0.40 


902 


43 


5.1±2.5 


0.19 


1 1.8 


22 


0.00 


0.00 


35.80zt0.67 


823 


232 


3.22±0.50 


839 


21 1 


1 l.I±1.9 


0.09 


71.6 


47 


0.49 


-0.1 1 


28.95zt0.52 


865 


204 


<2.52 






>11.5 




57.9 




0.98 


-0.22 


9.95zt0.77 


962 


150 


2.52±0.60 


986 


1 10 


3.9±1.2 


0.29 


19.9 


13 


1.96 


-0.44 


3.77zt0.92 


998 


88 


<1.76 






>2.1 




7.5 




2.94 


-0.66 


2.13zt0.64 


1030 


72 


< 1.18 






>1.8 




4.3 




-0.49 


0.1 1 


14.06zt0.46 


751 


177 


3.07±0.49 


745 


206 


4.6±0.9 


0.25 


28.1 


16 


-0.98 


0.22 


16.16ztl.47 


738 


1 13 


<2.04 






>7.9 




32.3 




-1.96 


0.44 


6.45zt0.81 


658 


130 


1.35±0.48 


675 


83 


4.8±2.3 


0.23 


12.9 


17 


-2.94 


0.66 


2.95zt 1.1 1 


636 


64 


<2.22 






>1.3 




5.9 




0.00 


0.00 


17.36zt0.62 


623 


144 


1.65±0.46 


606 


146 


10.5±3.3 


0.10 


34.7 


44 


0.51 


0.02 


23.77zt0.24 


674 


93 


2.53±0.46 


677 


80 


9.4±1.8 


0.1 1 


47.5 


39 


1.00 


0.03 


13.82zt0.67 


673 


152 


1.48±0.42 


697 


198 


9.3±3.0 


0.1 1 


27.6 


39 


1.50 


0.06 


4.98zt0.26 


726 


102 


0.71 ±0.27 


704 


108 


7.0±3.0 


0.15 


9.9 


28 


2.00 


0.06 


4.07zt0.25 


715 


1 12 


1.25±0.21 


731 


148 


3.3±0.7 


0.37 


8.1 


10 


2.50 


0.1 1 


3.29zt0.19 


718 


90 


<0.52 






>6.3 




6.6 




2.99 


0.22 


2.48zt0.33 


739 


70 


<0.80 






>3.1 




5.0 




-0.51 


-0.02 


19.50ztO.27 


571 


105 


1.76±0.31 


528 


141 


1 1.1±2.1 


0.09 


39.0 


47 


-1. 00 


-0.03 


1 1.40zt0.44 


590 


146 


1.77±0.47 


600 


140 


6.4±1.9 


0.17 


22.8 


25 


-1.50 


-0.06 


3.41zt0.23 


524 


99 


<0.83 






>4.0 




6.8 




-2.00 


-0.07 


1.89ztO. 16 


529 


108 


<0.52 






>3.6 




3.8 




-2.50 


-0. 1 1 


1.96zt0. 14 


535 


83 


<0.40 






>4.9 




3.9 




-2.99 


-0.22 


2.20zt0.31 


526 


76 


<0.90 






>2.4 




4.4 




-3.50 


-0.38 


2.86±0.19 


521 


51 


<1.10 






>2.6 




5.7 




-3.68 


-0.30 


3.68±0.30 


510 


43 


0.71 ±0.25 


486 


34 


5.2±2.2 


0.21 


7.4 


19 


0.00 


0.00 


8.01 ±0.70 


340 


153 


1.12±0.29 


391 


114 


7.2±2.5 


0.15 


16.0 


28 


0.00 


0.00 


16.15±1.05 


527 


214 


2.09±0.86 


527 


275 


7.7±3.7 


0.14 


32.3 


31 


0.00 


0.00 


33.35±0.39 


475 


103 


3.86±0.35 


472 


87 


8.6±0.9 


0.12 


66.7 


35 


0.00 


0.50 


28.90±0.42 


434 


80 


3.01±0.48 


425 


37 


9.6±1.7 


0.11 


57.8 


40 


0.00 


1. 00 


17.45±0.56 


410 


30 


2.43±0.41 


408 


23 


7.2±1.4 


0.15 


34.9 


29 


0.00 


1.50 


10.35±0.50 


406 


30 


1.26±0.41 


409 


41 


8.2±3.1 


0.13 


20.7 


34 


0.00 


2.00 


9.70±0.45 


401 


32 


1.98±0.59 


406 


26 


4.9±1.7 


0.23 


19.4 


18 


0.00 


2.50 


6.90±0.37 


398 


29 


0.80±0.29 


396 


37 


8.6±3.5 


0.12 


13.8 


36 


0.00 


3.00 


2.62±0.56 


388 


20 


<1.00 






>2.6 




5.2 




0.00 


-0.50 


27.35±0.53 


513 


62 


2.28±0.49 


520 


58 


12.0±2.8 


0.09 


54.7 


51 


0.00 


-1.00 


16.25±0.55 


534 


49 


1.90±0.44 


536 


42 


8.5±2.3 


0.12 


32.5 


35 


0.00 


-1.50 


6.55±0.23 


534 


51 


0.86±0.22 


563 


53 


7.6±2.2 


0.14 


13.1 


31 


0.00 


-2.00 


7.90±0.37 


538 


36 


1.64±0.38 


539 


43 


4.8±1.3 


0.23 


15.8 


17 


0.00 


-2.50 


5.77±0.33 


534 


46 


<0.68 






>8.5 




11.5 




0.00 


-3.00 


2.28±0.35 


553 


36 


<0.62 






>3.7 




4.6 








11.85±0.30 


462 


151 


1.73±0.37 


411 


20 


6.8±1.6 


0.16 


23.7 


27 


0.00 


0.00 


9.82±0.51 


255 


72 


1.97±0.49 


277 


128 


5.0±1.5 


0.22 


19.6 


18 



NGC 2903 
NGC 3031 
M82 



NGC 3521 
NGC 3627 



NGC 4736 
NGC 5055 
M51 



M51 disk" 
NGC 5457 



° Offset from the nucleus position listed in Table[T] in units of ai'cminutes. 

' The measured integrated intensities and associated uncertainties, calculated using the prescription explained in the text. For 
non-detections, a 2 cr/ upper limit was given. 

Velocity and line widths are Gaussian fit values, or else are calculated from the moment for non-Gaussian lines. 

The ratio of ^^CO to ^'^CO integrated intensities. The errors are based on the statistical uncertainties of integrated intensities, 

which can be derived from the error transfer formula o"(7^i2/i3) = ([^^~^]'^ + [ ^ ^'^^Ji3 ^^^^^ j^)-*-/-^. 
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Table 3 Observed and Derived Properties of 



Source 


Aa 


AS 


/clSo ± cri 


Vmcan 


AV 




N{H2 




Source 


arc mi n 


arcmin 


(K km s^^) 


(km s ) 


(km s ) 






M82 


0.0 


0.0 


3.81±0.90 


140 


121 


0.017 




18.6 




0.0 


0.5 


<4.2 














0.0 


-0.5 


2.55±1.27 


224 


176 


0.022 




14.0 




0.5 


0.0 


3.80±0.90 


145 


93 


0.011 




25.2 




0.5 


0.5 


3.30±0.76 


289 


135 


0.012 




20.5 




0.5 


-0.5 


<4.1 














1.0 


0.0 


4.19±1.20 


284 


236 


0.018 




20.6 




1.0 


0.5 


3. Mil. 84 


256 


174 


0.015 




14.8 


M51 


0.0 


0.0 


1.04±0.35 


466 


176 


0.032 




7.4 




0.0 


0.5 


0.92±0.49 


413 


154 


0.032 




5.7 




0.0 


1.0 


0.45±0.16 


407 


23 


0.026 




4.6 




0.0 


1.5 


0.99±0.54 


393 


79 


0.10 




2.4 




0.0 


2.0 


<1.17 














0.0 


2.5 


0.88±0.37 


369 


58 


0.14 




1.6 




0.0 


-0.5 


0.95±0.17 


519 


37 


0.035 




4.3 




0.0 


-1.0 


0.55±0.13 


544 


35 


0.034 




3.6 




0.0 


-1.5 


<0.68 














0.0 


-2.0 


<0.84 













The average optical depth in C^°0 emission hne is calculated from the similar equation as eq.(4). 
H2 column density derived from C^*0, calculated from eq.(6). 
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Table 4 Line Intensity Ratios 



Source 


Aa 
arcmin 


AS 
arcmin 


-fi^coZ-'^ciSo 


■fi3co/-^ci8o 


^13co/^HCO + 




References 


M82 


0.0 


0.0 


53.5±13.0 


2.7±0.9 


1.4±0.2 


2.1 ±0.5 


This work 




0.0 


0.5 


>35.3 


>2.8 


0.9±0.1 




This work 




0.0 


-0.5 


40.5±20.6 


2.2±1.4 


0.5±0.1 




This work 




0.5 


0.0 


82.3±19.5 


3.7±1.2 


1.2±0.2 


3.0±0.9 


This work 




0.5 


0.5 


76.5±18.0 


3.4±1.1 


0.9±0.1 


4.0±1.4 


This work 




0.5 


-0.5 


>29.8 


>1.4 


1.6±0.7 




This work 




1.0 


0.0 


49.5±14.5 


2.7±1.1 


1.8±0.4 


3.3±1.0 


This work 




1.0 


0.5 


63.2±37.5 


2.0±1.4 


1.9±0.4 




This work 


NGC 3628 


0.0 


0.0 






0.5±0.1 




This work 


NGC4631 


0.0 


0.0 


>22.8 


>2.4 






This work 


M51 


0.0 


0.0 


29.1±10.3 


3.7±1.6 


0.9±0.2 


2.6±0.9 


This work 




0.0 


0.5 


28.6±15.8 


3.3±2.3 


1.8±0.7 


2.4±0.8 


This work 




0.0 


1.0 


34.9±13.7 


5.3±2.8 


1.3±0.7 




This work 




0.0 


1.5 


9.4±5.6 


1.3±1.1 






This work 




0.0 


2.0 


>7.6 


>1.7 






This work 




0.0 


2.5 


7.1±3.4 


0.9±0.7 






This work 




0.0 


-0.5 


26.0±5.2 


2.4±0.9 


1.5±0.5 


5.1±2.1 


This work 




0.0 


-1.0 


27.0±7.1 


3.5±1.6 






This work 




0.0 


-1.5 


>8.7 


>1.3 






This work 




0.0 


-2.0 


>8.6 


>2.0 






This work 


M51 disk 






16.4±7.4 


2.6±1.7 






This work 


NGC 253 


0.0 


0.0 


61.0±6.2 


4.9±0.7 


1.8±0.1 


2.5±0.4 


1.2,3 


NGC 891 


0.0 


0.0 








4.5±1.4 


3,4 


NGC 1068 


0.0 


0.0 








1.7±0.5 


3,5 


NGC 1808 


0.0 


0.0 


48.9±5.4 


3.3±0.6 






6,7 


NGC 2146 


0.0 


0.0 


32.7±9.1 


3.0±1.1 






6,7 


NGC 2903 


0.0 


0.0 








4.5±1.5 


3,4 


NGC 3256 


0.0 


0.0 


91.6±28.1 


2.9±1.5 


0.8±0.3 


1.8±1.1 


8 


NGC 4736 


0.0 


0.0 








4.5±1.5 


3,4 


NGC 4826 


0.0 


0.0 


19.0±2.7 


4.1±1.1 






6,7 


NGC 5457 


0.0 


0.0 








2.9±0.8 


3,4 


NGC 6240 


0.0 


0.0 






0.2±0.1 




9 


Maffei2 


0.0 


0.0 








5.2±0.8 


3,4 


IC 342 


0.0 


0.0 


35.1±3.0 


3.9±0.3 




2.1 ±0.4 


1,3,4 


Circinus 


0.0 


0.0 


54.4±5.4 


5.7±0.7 






6,7 


Aip 220 


0.0 


0.0 


49.4±19.9 


1.0±0.4 


0.3±0.1 


0.6±0.3 


9 



Notes: All the rati os of integra ted intensity have been corrected for the different beamsizes. 

References: (11 ISage et al. I il99lh; (2) iHenkel et al. I il993h: (3 ) ISage et al. I il990l): (4) ISage & Isbell il99lh : (5) 
lYoung & Sandersi i ll986l) : f6) lAalto et al1il99ll) : (7l lAalto et al. Hl995l) : fSl lCasoh et al. Ml99a : t9l lGreve et al. U2009I) . 



26 



Tan et al. 



Appendix A: THE STABILITY OF DLH RADIO TELESCOPE 



Figure lATTI is an Allan Variance Plot, which is often the ultimate way to measure stability but requires 
an enormous amount of observation time. There are three main contributions to be aware of in the 
Allan plot, including the white noise, the l//-noise and low frequency drift noise. The upper panel of 
Fig. lA. II shows that how the squared RMS noise of ^^CO spectra with velocity resolution of ~0.16km 
s^^ varied with the integration time on source, while the lower panel shows the relative error deviate 
from the radiometer equation. The radiometer equation, also the limiting sensitivity of the spectrometer, 
is given by 

^ = ^, (Al) 

Tsys V AVT 

where Tgyg is system temperature, r is the sum of the integration time on the source and on the off 
position, K is a factor accounting for data taking procedures, and ATmis is rms noise temperature of 
the spectra for a given frequency resolution Ai/. 

Figure lA.lfa and b show the data obtained from the observing semesters of 2008 to 2009 and of 
2009 to 2010, respectively. It can be seen in Figure lATh that the Allan plot begin to deviate from the 
radiometer equation when integrating about 10 minutes on source, and the relative error increase to 80% 
when integrating 100 minutes, however, it is shown in Figure lATb that the relative error increase to only 
10% with the same integration time. Therefore, both the stability and the sensitivity of the telescope have 
been greatly improved after the system maintenance in the summer of 2009. The average noise level in 
units of r| would be 0.020 K and 0.018 K when integrating 100 and 200 minutes on source respectively, 
with spectra velocity resolution of ~^10km s^^. Therefore, this plot can be used to estimate observation 
time according to the sensitivity that we required. But an important thing to note in this Allan Variance 
Plot is that the data we used are the raw data, that is to say, we didn't do any data processing such as 
rejecting the spectra with distorted baseline or abnormal rms noise level. Consequently, we could get 
even better sensitivity if we just co-added the spectra that found to be consistent within normal rms noise 
level. Based on these analysis, it seems to imply that the telescope still have the capability to detect even 
weaker emission signal than what have been detected, since the effective integration time of most of our 
observations have not reach the limit. 
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Fig. A.l The Allan Variance and the relative error from radiometer equation as a function of 
integration time on source (see Appendix |A}. On the upper panel of each figure, the dashed 
line represents what is expected from the radiometer equation with a slope of -1. On the lower 
panel, each point represent the relative error between the value of measured and expected 
from radiometer equation, (a)the data taken from the observing semester of 2008 to 2009. 
(b)the data taken from the observing semester of 2009 to 2010. 



